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Abstract 



Using methods developed for hard exclusive QCD processes, we calculate 
the asymptotic behaviour of heavy-meson form factors at large recoil. It is 
determined by the leading- and subleading- twist meson wave functions. For 
1 \v ■ v'\ <^ mq/A, the form factors are dominated by the Isgur-Wise 
function, which is determined by the interference between the wave functions 
of leading and subleading twist. At |w • f'| ^ mq/A, they are dominated 
by two functions arising at order 1/mq in the heavy-quark expansion, which 
are determined by the leading-twist wave function alone. The sum of these 
contributions describes the form factors in the whole region |f • t;'! ^ 1. As 
a consequence, there is an exact zero in the form factor for the scattering 
of longitudinally polarized B* mesons at some value v ■ v' ^ mi,/ A, and an 
approximate zero in the form factor of B mesons in the timelike region [v-v' ~ 
—rrib/A). We obtain the evolution equations and sum rules for the wave 
functions of leading and subleading twist as well as for their moments. We 
briefly discuss applications to heavy-meson pair production in e~^e~ collisions. 

12.38.Bx, 12.39.Hg, 13.40.Gp, 14.40.Nd 
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Typeset using REVT^ 



I. INTRODUCTION 



In the heavy-quark effective theory (HQET) (see for reviews), the heavy- 

quark spin does not interact with gluons to leading order in 1/ uiq (where mg is the heavy- 
quark mass). Therefore, this spin can be rotated (spin symmetry) or even switched off 
(superflavour symmetry fT^JTSU ) without affecting the dynamics. In the heavy-quark limit, 



the properties of the doublet of the ground-state pseudoscalar and vector mesons (qQ) are 
therefore characterized by the spin-parity quantum numbers = of the light degrees 
of freedom |l;14]. In this paper, we shall use the superflavour symmetry to describe the 



ground-state mesons by a Dirac wave function. However, we collect in Appendix A the most 
important formulae using a more conventional formalism. 

Let Ql be a scalar field describing a heavy antiquark moving at four-velocity v with 
its spin switched off. Then the decay constant / of a heavy meson (moving at the same 
velocity) is defined as 

{0\Q:q\M{v)) = fu{v), (1.1) 

where u{v) is the Dirac wave function of the meson, which satisfies 

'^u{v) = u{v). (1.2) 

A non-relativistic (i.e. mass independent) normalization of u{v) and \M{v)) is assumed. In 
the heavy-quark limit, the relation between / and the usual meson decay constants reads 

/a/ = /m- = ^. (1.3) 

To leading order in l/mg, current-induced transitions between two ground-state mesons 
are described by a single Isgur-Wise form factor |lll,|l5| : 

{Miv')\Q:Q,,\Miv))=avv')uiv')uiv), (1.4) 

where v, v' are the meson velocities. At next-to-leading order, there appear l/mq corrections 
to the currents and to the Lagrangian of the HQET |]16[. The first type of corrections can 



be expressed via the matrix element of a dimension-four operator:^ 

{M{v')\{iD^^Q:)Q,,\M{v)) = u{v')e{v,v')u{v), (1.5) 

where 

^^'{v, v') = ^+{v ■ v'){v + v'Y + ^-{v- v'){v - v'Y + ^siv ■ v') Y ■ (1-6) 



The equations of motion, iv ■ DQ^ = 0, can be employed to relate ^± to ,^3 [0. The result 
is: 



^We use Di^ = - iA^" and D^t = 5^ + iA^^\ where A^" = QstaAi^ is the gluon field. 
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(a) (b) 

FIG. 1. Hard-gluon exchange contributions to heavy-meson form factors. The external current 
is presented by the wave line; the heavy antiquark is represented by a double line. 



V ■ v' + 1 



[1.7] 



where A is the "binding energy", i.e. the difference between the meson mass and the heavy- 
quark mass. Hence, there is only one new independent form factor. The l/rriQ corrections 
to the Lagrangian give rise to the matrix elements of two non-local operators: 



{M{v')\ I j dxT{Q:Q,,(0), Ql{iDfQ,{x)}\M{v)) = 2xi{v ■ v') u{v')u{v) , 

{M{v')\i I dxT{QlQ,\{)),QliG^''Q,{x)}\M{v)) = 2u{v')x''''{v.v')u{v) , (1.8) 



where iG^"" = [iD^,iD''] = igstaG^" is the gluon field strength, and (af"" = 

X'^'^iv, v') = X2{v ■ v'){Yv"' - Yv'n + 2tX3{v ■ vy^" . (1.9) 

In the second matrix element in (|1.8|), the indices /x and z/ are restricted to the subspace 
orthogonal to v. Explicit expressions for the meson form factors in terms of the functions 
^, ^3 and Xi can be found in |l9|,[l6[ 



At moderate values of v 



these functions have 



been studied extensively in the framework of QCD sum rules, both at leading [p!^^ -|2T|] and 
next-to-leading |2^-p4| order in the heavy-quark expansion. In the present paper, we shall 
consider the behaviour of the form factors in the large-recoil region |f ■ f '| ^ 1. This region 
is inaccessible in the weak semileptonic decays, but it can be explored (at least in principle) 
in the production reaction e"'"e~ — ^ M^*^M^*\ Using methods developed for hard exclusive 
processes, we calculate the asymptotic behaviour of the form factors in a model-independent 
way. 

Our results can be summarized as follows: For |f ■ f'| ^ 1, there is a large momentum 



transfer to the light quark: 



ight 



-A^ 



V ■ V 



where A is of the order of a typical hadronic 



mass scale. As shown in Fig. |T]a, this momentum is transferred by the exchange of a hard 
gluon, and the methods developed for hard exclusive processes in QCD [^5H 28|| (see [^^ -|3l[ 
for reviews) are applicable. In the "brick wall" frame, where v' = —v, the projection of the 
total angular momentum on the z axis (directed along v) is equal to the projection of the 
meson spin. (Recall that the heavy-quark spin has been switched off.) Since it is conserved, 
the meson helicity changes its sign |^|. The asymptotic behaviour of the Isgur-Wise form 
factor is thus determined by the interference between the leading-twist (quark-antiquark) 
wave function and the subleading-twist (quark-antiquark and quark-antiquark-gluon) wave 
functions. It is given by 
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^(^ ■ ~ TI^^ ■ (1-10) 

Indeed, the situation is similar to the well-known case of the vr-p form factor pOf. At 
order l/mg in the heavy-quark expansion, there are however contributions involving the 
leading-twist wave function only. They conserve the meson helicity and behave as 

^^{v-v') v-v'x2{v-v') asp , . 

~ ~ . I i. ii I 

niQ tuq itlqA'^ V ■ v' 

The leading contribution in the heavy-quark expansion, which is given by the Isgur-Wise 
function, dominates as long as \v ■ v'\ <^ mg/A. For \v ■ v'\ ^ rriQ/A, however, the con- 
tributions of ^3 and X2 in ( |1 . 1 1| ) become the dominant ones. Note that they violate the 



heavy-quark spin symmetry, i.e. they contribute in a non-universal way to the various meson 
form factors. Higher-order terms in the heavy-quark expansion (of order l/mg and higher) 
cannot fall off slower than l/(w -v') because this behaviour corresponds to the leading twist, 
and hence they always remain small corrections. 

It is instructive to consider the same situation from an opposite point of view. At 
asymptotically large values ~ 2mQ\v-v'\ ^ itt-q/x (where q is the momentum transferred 
to the mesons, and x ~ A/mg is the momentum fraction carried by the light quark), the 



form factor of a heavy meson behaves like that of the pion [25-31]: 



F{<l')-^, (1.12) 



which exactly corresponds to ( |1 . 1 1|) . However, there is a contribution to the subleading-twist 
(1/g^) correction which is proportional to rriQ. It becomes important for moderate values of 



F(g^)-^^. (1.13) 

This contribution exactly corresponds to ( p..l(]| ). It dominates for jg^l -C rnqX. Higher- twist 
(1/g^ and higher) corrections cannot be more enhanced than tuq because otherwise the 
form factor would diverge in the heavy-quark limit, and hence they always remain small 
corrections. 

Until now, we considered form factors for transitions induced by a current containing 
heavy quarks only. In the case of, say, the electromagnetic current, which also contains 
light-quark fields, contributions of the type shown in Fig. |l]b appear. However, they lead to 
the behaviour 

- ^ (1.14) 

UIq V ■ v' 



and can thus be safely neglected, since rriQ ^ A for a heavy quark. 

In summary, for 1 <^ |f ■ f'| ^ mq/ K the dominant contribution to meson form factors 
comes from the universal Isgur-Wise function. It involves the exchange of a longitudinal 
hard gluon and a change of the meson helicity, corresponding to the interference between the 
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leading- and subleading-twist wave functions. For |f ■ f'| ^ itiq/A, on the other hand, the 
situation is the same as in massless QCD. The asymptotic behaviour of meson form factors 
is determined by leading-twist contributions, which are not universal. They are governed 
by the exchange of a hard transverse gluon, which conserves the meson helicity. The sum 
of these two contributions describes the asymptotic behaviour of the form factors in the 
whole region |f ■ f'| ^ 1. The simple picture described here is only slightly modified by 
the emission of gluon bremsstrahlung, which can be dealt with in renormalization-group 
improved perturbation theory. It leads to an additional, moderate power suppression of the 
form factors at large recoil. 

The remainder of the paper is organized as follows: In Sec. we introduce the (quark- 
antiquark) meson wave functions of leading and subleading twist, ip±{uj). Unlike in QCD, 
wave functions defined in the HQET depend on a dimensional argument u. We investigate 
the moments of these wave functions and derive the symmetry relations between the vari- 
ous meson wave functions, which arise in the heavy-quark limit. In Sec. |T|, we derive the 
evolution equations for ip±{u), which are analogous to the Brodsky-Lepage equations p8| 



This calculation extends the calculation of the HQET anomalous dimension of local heavy- 



light current operators [p^ -|38[| . A new kind of ultraviolet divergence appears in the relation 
between the local operators and the operators defining the wave function. Therefore, the 
Brodsky-Lepage kernels do not determine the renormalization properties of the local op- 
erators completely. A similar situation is encountered in the case of the (Altarelli-Parisi) 
equations describing the evolution of distribution functions in the HQET ||39|. In Sec. we 
investigate the properties of the wave functions y:>±{uj) using the QCD sum-rule approach. 
This extends the heavy-meson sum rules |p!8| , ^0| -^ to the case of non-local operators. After 



considering the sum rules for the lowest moments, we construct the sum rules for the wave 
functions themselves, taking into account the non-locality of the quark condensate [Q. In 
Sec. 0, we apply our results to derive the leading asymptotic behaviour of meson form 
factors at large recoil. First, we calculate the contribution of the interference between the 
leading-twist and the (quark-antiquark) subleading-twist wave functions to the asymptotic 
behaviour of the Isgur-Wise function. | Then we calculate the leading-twist contributions 
to the form factors appearing at order l/mq in the heavy-quark expansion. Finally, in 
Sec. VI we discuss the implications of our results for the reactions e+e- ^ and 



e+e" Technical details of our calculations are presented in four appendices. 

Before we proceed, some comments on the existing literature on the application of pertur- 
bative QCD to the calculation of heavy-meson form factors are in order. A simple model of 
a meson made out of two heavy quarks with unequal masses was considered in |^^. There 



it was noted that some form factors must have zeros in the physical region. We confirm 
this interesting observation, although we disagree with some other results of this work (see 
Sec. 0). A similar model was considered in |]45|- ^] . There, a single spin structure of the 



heavy-meson wave function was used, which is determined from the condition that the light 
quark be at rest in the meson rest frame. Hence, all quark-antiquark wave functions have 
the same shape, and this shape was assumed to be 6{lj — fi), with /i being the constituent 



^The properties of quark-antiquark-gluon subleading-twist wave functions and their contribution 
to ^,{v ■ v') will be discussed elsewhere. 
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mass of the light quark. As we shall see later, for a realistic heavy meson the wave functions 
(p+{uj) and (^-(cu) do not coincide, and they are not well approximated by sharply peaked 
functions. Integrals for the form factors receive important contributions from the region 
of low uj values, which are missing in the peaking approximation. Therefore, the results 
obtained using such a static quark model can at best be taken as a crude estimate. Pertur- 
bative QCD and the constituent quark model were recently applied also to the semileptonic 
decays B — > D^*H u H3], for which 1 < f ■ f' < 1.6. In our opinion, such small values of v -v' 



are far too low to treat the gluon with ~ —A v ■ v' perturbatively. Moreover, the calcu- 
lations in pSf are done using a model wave function with an ad hoc dependence, whose 



longitudinal-momentum dependence contradicts the expectations based on the HQET. 

II. QUARK ANTIQUARK WAVE FUNCTIONS 

We shall define the quark-ant iquark wave functions ^±{t) of a heavy meson in terms of 
the matrix element of the bilocal operator 

d{t)=Q*{<<))E{<d,z)q{z)- t = vz, (2.1) 

where 2; is a null vector on the light cone {z'^ = 0), and 

E{x,y)=Pexp(^-iJdz^A^{z)^ (2.2) 

is a string operator ensuring gauge invariance. In the light-cone gauge (^4+ = 0), one simply 
has E{0, z) = 1. Since in this section we are considering operators containing a single heavy 
quark field Q*, we shall for simplicity omit the velocity label on the field. Similarly, we shall 
write M and u instead of M{v) and u{v). The meson matrix element of the operator 0{t) 
has two independent Dirac structures, u and ^ u, and we define 

( |0(t)|M) = / (^+(t) + 1 [^4t) - ^+(t)] ^) u . (2.3) 

It is convenient to introduce two light-cone vectors ra^ = (1,0,0,=f1) such that = 
and ■ n_ = 2. Any vector can be decomposed as a'^ = ^(a^n^ + a^nl) + aj_, where 
a± = a ■ n±. This implies a ■ b = |(a+6_ + a_6+) — a± ■ b±. We shall also use the light-cone 
components of the Dirac matrices, defined as 7-1- = n±7^ = 1^±- If the meson is at rest, then 
v'^ = |(n+ + n'^), i.e. = v_ = 1. Using u = u, we can then rewrite ( p.3|) as 

(0 |d(t)|M) = i/ [^+(t)7- + ^-ith+\u. (2.4) 
For a meson with an arbitrary velocity in the plane, this formula becomes 

{0\d{t)\M) = \f [^+(t)i;+7- + <?-(i)^-7+] (2.5) 
If we introduce the rapidity 'd by writing 

w'' = (cosh ■(9,0,0, sinh 7?) , (2.6) 



6 



then v+ = and f_ = e~^. This shows that for a fast-moving meson (-i? ^ 0), is the 
leading-twist wave function, whereas has subleading twist. It is convenient to project 
onto these wave functions by writing 

( |0±(t)|M) = / ^±(t) 7±M , 0±(t) = 7±0(t) . (2.7) 

This result is valid in an arbitrary reference frame, as can be seen by using the relations 

7i = 0, 7±7t = — 7±/'- (2.8) 

v± 

The wave functions (p±{t) depend on the separation t on the light cone. We define the 
corresponding wave functions in momentum space by 

<^±M = ^J dt^i{t)e''^\ <^±(t) = J duM^)e-''^'. (2.9) 

The variable u has the meaning of the light-cone projection of the light-quark momentum 
in the heavy-meson rest frame. The positions of the singularities in the complex t plane are 
such that ^p±{u!) vanish for uj < 0. The wave functions are normalized such that 

oo 

(^(0) = Jdujip±{uj) = 1. (2.10) 



We can formally introduce operators 0±{uj) such that 

(0|O±(a;)|M) = /(^±(a;)7±M. (2.11) 

This implies 

0±i^) = ^J d'^^±(t) e^-* = g*(0)7±5(^/^+-u;)g(0), 
0±it) = J duj 0±{uj) e-'""' . (2.12) 

Expanding in powers of t in the definitions ( p.9|) and ( p.l2|) , we obtain 

o±(t) = E or ^ , ^±it) = E , (2.13) 



„=o "o 



where 



Ol") = J dujO±{uj)uj" = Q*l±{tD+Yq, 
(w")± = / dw(^±(cj)cj". (2.14) 



Equation ( p.ll ) then implies a relation between the moments of the momentum-space wave 
functions and the local, higher- dimensional operators O^'^: 



0|Oi")|M) = /(cu")±7±«. (2.15) 
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Using the equations of motion, the first moments of the wave functions can be calculated 
in terms of the parameter A encountered in (|1.7| ) In general, we may write 

(2.16) 



; \Q*iD^q\M) = f {av^" + b-f'')u . 



The equations of motion for the light quark, ipq = 0, imply that (a + 46) = 0. The equations 
of motion for the heavy quark, iv ■ DQ = 0, can be used to write 



( \Q*iv ■ D q\M) =ivd{Q \Q*q\M) = A ( |g*g|M) 



(2.17) 



where A = rriM — rriq is the effective mass of the meson M in the HQET ||50|. This relation 
implies that {a + b) = A, and therefore 



( \Q*iD''q\M) = i/A {Av^" - 7^)m , 



3,..,,. , (2.18) 
Using this result, we find that the first moments of the wave functions are given by 



^A. 



(2.19) 



A similar analysis can be performed for the second moments. Consider the matrix 
element 



{0\Q*iDf'iD''q\M) = fe^^u, 
where the most general form of 6^^^ is 



(2.20) 



(2.21) 



The equations of motion impose three independent relations among the five parameters q, 
which imply that the matrix element in (|2.20| ) is completely determined by its antisymmetric 
part i.e. by the matrix element of the gluon field iG^'^ = [iD^,iD'^]. For reasons to 
become clear below, we find it convenient to introduce two hadronic parameters A|; and 
by 



C5 



1 \2 



H ■ 



^4 — |(A^ - X%] 
In terms of these quantities, we obtain 

( mC^'qlM) = \f - XDiYv" - Yv'') + iXla^-" 



(2.22) 



u . 



: |Q*H^^^ ^D-"} q\M) = \f [(6A2 + 2A| + \l)v^v'' - (A^ + A| + \^h)9 



u . 



(2.23) 



From the second relation, it follows that the second moments of the wave functions are given 
by 



2A^ + |A| + ^X^ , 



1 \2 



2X2,1x2 
3A + ^Ajj . 



(2.24) 



According to the first equation in (|2.23| ), the moments {uj'^)± are thus related to normaliza- 
tion integrals of quark-antiquark-gluon wave functions. 
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Our definition in ( |2.22| ) is such that, in the rest frame of the heavy meson, the quantities 
A|; and parametrize the matrix elements of the chromo-electric and chromo-magnetic 
fields, respectively. Defining^ Ei = Goi, Hi = —^eijkGjk, and = 7^7*, we find 



{0\Q*ia- Eq\M) = fX^u, -{0\Q*a ■ H q\M) = fX^u. 



(2.25) 



To finish this section, let us switch the heavy-quark spin on and relate the numerous 
quark-ant iquark wave functions of the ground-state pseudoscalar and vector mesons, M 
and M*, to the HQET wave functions ip±. These relations are very conveniently obtained 
using the covariant tensor formalism described in Appendix A. For a pseudoscalar meson 
M, the matrix elements of the pseudoscalar, axial, and tensor currents are non-zero, and we 
define a set of four wave functions in the following way:0 



(0|Q(0)75gW|M) = -ifMmM^p, 
{ |Q(0) 7^75 qiz)\M) = fM^AiP" - rriM ^A2 
: |Q(0) ^^^^75 g(^) |M) = ifu {P^z^ - P'z^) , 



(2.26) 



where (pi = ipi{p- z) = (piimMt)- For simplicity, we have omitted the string operator £'(0, z), 
i.e. we have adopted the light-cone gauge = 0. In the heavy-quark limit, we obtain, 
using the results of Appendix A: 



VP 



VAX = V+if) 5 VA2 = VT 



2 



t 



For t = 0, we obtain in this limit the normalization conditions: 

(?P(0) = (^Al(0) = 1, VA2{0) 



<^t(0) = - 



(2.27) 



(2.28) 



where the second relation is a consequence of (|2.13|) and ( p.l9|) . 

For a vector meson M* with polarization vector e, the matrix elements of the scalar, 
vector, axial, and tensor currents are non-zero, and we introduce a set of six wave functions 
as follows: 



( |g(0) g(^)|M*) = ifM-mu' Vsz-e, 
{ |g(0) qiz)\M*) = Jm* \mM' Vvi + IVV2 {p-ze^ 
(0|Q(0)7^5g(^)|M*) 
(0|Q(0)a^^g(;^)|M*)=/M 



e ■ zp^ 



i^Ti {e^p" - eV) 



rriM* VT2 {e^z'' - z^) 



(2.29) 



^If we define = (D^, -D), then E = i[D^,D] and H = -iD x D. 

'^Contrary to the notation used in the rest of this paper, here and in ( p. 29 ) we use the stan- 
dard relativistic normalization of states, which adds a factor ^mu on the right-hand side of the 
equations. 
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In the heavy-quark hmit, we find the relations 



~ / \ ~ ~ ~ ~ 'i' — ^-(t) 

^Vl = VTI = 'P+yt) , = fV2 = fA = fT2 = - 



and the corresponding normahzation conditions 

m(o) = <^Ti(o) = 1 



A 



(fsiO) = V?V2(0) = (fAiO) = V?T2(0) = - 



The QCD wave functions in momentum space are defined as 

1 



- /d(p■z)^.(p■^)e-"•^■^ 



(2.30) 



f2.311 



(2.32) 



so that X = uj/rriM- On the basis of the behaviour of the eigenfunctions of the evolution 
equations as well as sum-rule inspired arguments, it is usually assumed that for pseudoscalar 
mesons ipAi{x) ~ x and (pp{x) ~ 1 as x — For the HQET wave functions, this implies 
the behaviour 



(2.33) 



as — > 0. In Sec. |V|, we will indeed find these scaling laws from an explicit calculation of 
the wave functions using QCD sum rules. 



III. EVOLUTION EQUATIONS 

The definitions of the previous section are somewhat formal, because the operators in- 
volved require renormalization. In this section we discuss how the ultraviolet divergences 
in operator matrix elements can be removed in a consistent way. After doing this, how- 
ever, we shall ignore renormalization effects in the further course of this paper. The reader 
not interested in the conceptual problem of renormalization can thus proceed directly with 
Sec. 0. 

We use the MS subtraction scheme in d = 4 — 2e dimensional space-time. The bare and 
renormalized operators are related by 

Ol^'^Lu) = J duj' Z±{uj, to') 0±(cu') , (3.1) 

where 

Z± (cu, u') = 6{uj -u')- — z± {u, uj') + ... , (3.2) 

and the ellipses represent poles of higher order in 1/e. The operators 0±{uj) and hence their 
matrix elements fip± (uj) obey the renormalization-group equations 



dlog /i 

where the anomalous dimensions r±{uj,uj') are given by 



J did' T±{uj, uj') f(p±iuj') = , (3.3) 
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(b) (c) 

FIG. 2. One-loop diagrams contributing to the matrix elements {0\0^'^'^{oj)\ll!'). The bare 
current operators are represented by a circle. 

d 

T±{uj,uj') = asT:— z±{uj,uj') . (3.4) 
da. 



Equation ( |3.3| ) is analogous to the Brodsky-Lepage evolution equation in QCD 

In order to obtain the anomalous dimensions T±{uj, uj') at one-loop order, we consider the 
matrix elements (0 |0^'^'''^(ci;)|a;'), where \uj') represents a state consisting of a scalar heavy 
antiquark at rest and a light quark with momentum = u'. According to (|3.1| ), these 
matrix elements equal Z±{u!,uj') 'y±u. The relevant one-particle irreducible loop diagrams 
are shown in Fig. ^ Although the operators 0±{uj) in ( p.l2|) take a particularly simple 
form in the light-cone gauge ^4+ = 0, we refrain from adopting such a singular gauge and 
work instead in the Feynman gauge.0 To obtain the Feynman rules for vertices involving the 
operators 0±{uj), we start from their definition as the Fourier transforms of the non-local 
operators 0±{t) and obtain, to first order in the gauge field: 

0±M = ^J dte^-*|Q*(0)7±e*^+g(0) 

-I ydrg*(0)7±e"^+A+(0)e(*-")^+g(0) + ...|. (3.5) 

It is then straightforward to derive the Feynman rules shown in Fig. ^. 

Let us then sketch the calculation of the diagrams in Fig. ^. The contribution of the first 
diagram is: 

zG^2^/i J clfc+dfc_ ^k^^^O)ik' + ^0)[{k-py + ^0] ' ^^'^^ 

where Cp = {N^ — l)/{2Nc). The integral is trivial to perform, and the k_ integral 
is readily calculated by the method of residues. The poles of the integrand are located at 
k- = —oj + iO, /c_ = {k^ — iO)/u!, and A;_ = {kj_ — iO)/{u! — u'). If c<j < 0, all poles lie in 
the upper half plane, and the integral vanishes. For c<j > 0, it is necessary to distinguish 
the cases u > u' and uj < oj' . For uj > cu', we close the contour in the upper half plane and 
set k_ = —UJ. For uj < uj' , we close the contour in the lower half plane and set /c_ = k\/uj. 
Only in the second case and for the minus projection there is an ultraviolet divergence, 



^We have repeated the calculation of Tj^{uj,uj') in the light-cone gauge and obtained the same 
result as in the Feynman gauge; however, we could not find an easy way to recover the correct 
result for r_(a;,u;'). 
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.A. . 




+ -P+ 



as ta n+ 7± ^ _ 



FIG. 3. Feynman rules for vertices involving 0±{lv) 



which arises from the fc_L integration. Keeping only the singular term, we obtain for the 
contribution to Z±{uj,uj'): 



(3.7) 



The other two diagrams in Fig. ^ are evaluated in a similar way. For the second one, we find 
an ultraviolet divergence for the plus projection if fc+ < uj' . Its contribution to Z±{uj,uj') is: 



a 



as 
Aire 



(1±1) 



1 



u' (id — u;') + 



9{lu' — lu) — 5{uj' — LU)> . 



The distribution l/{u — u')+ is defined such that 



duj 



duj 



u — u' 



(3.8) 



(3.9) 



for any smooth function f{uj). The third diagram has an ultraviolet divergence if > u'. 
Its contribution is: 



as 



1 



27fe (to — UJ 



(3.10) 



Finally, we have to add the contributions from the wave-function renormalization of the 
external lines, which gives [^] 



Z'^/'zy' 5{u - J) = (l + ^) biuj - J) . 

Collecting all terms, we obtain for the quantities z±{uj,u') defined in 

{lj' - uj) 9{uj - Lj') 3 



(3.11) 



z±{uj,u') =Cf — \ ± 

TT 



OJ' 



5{uj — oj') 



— 0{UJ — UO ) 



(3.12) 
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To one-loop order, the anomalous dimensions are given by the same expression. We thus 
obtain 



r+(uj,uj') = Cf 



a. 



TT 



1 ^ 0{lo' - LU) 



TT 



\uj — 1+ UJ' 

e{uj-Lu') e{iu'-uj) 



6{iU — Lu') 



{uj - UJ')+ 



UJ' 



— -6(uj — uj' 
4 



(3.13) 



In the familiar case of QCD wave functions for light-quark systems, the renormalization 
of the non-local operators analogous to 0±{u!) would suffice to renormalize the tower of local 



operators defined analogous to 0±^ in (|2.14| ). Accordingly, the renormalization of the wave 



functions v?i(x) renders their moments (a;")i finite. Unfortunately, this situation does not 
hold in the case of wave functions for heavy mesons defined in the HQET. The reason for 
this unexpected fact is that in the HQET the wave functions depend on the dimensional 
variable u, which takes values between and oo. As a consequence, equation ( |2.14| ), which 
defines the operators o£^^ in terms of weighted integrals of 0±{uj), does not hold for the 
renormalized operators; the integral over u has an additional ultraviolet divergence not yet 
removed by the renormalization of 0±(u;). This divergence must be removed separately. 

Consider, as an example, the simplest case n = 0. Then the bare operators Q^^'^^"^*^ = 
Q*1±Q are local heavy-light current operators with dimension three, which are renormalized 
in the following way: 



0(0), bare 



where 



To 



3 „ 

LjW — 

4 TT 



of 



(3.14) 



(3.15) 



is the well-known one-loop hybrid anomalous dimension [P5|-p5|. On the other hand, the 



bare operators can be expressed in terms of integrals over the renormalized operators 0±{uj) 
using (|0|). This gives: 



Q(0),bare 



j dudu' Z±{uj,uj')0±{uj') = (^j dcjZ±(cu,cu')) j dcuO±(cu) 



(3.16) 



As mentioned above, the last integral over the renormalized operators 0±{(jj) is ultraviolet 
divergent from the region uj oo, and we define 



j dujO^{uj) = z'^of = (^1 - Iri + . . .) o5 

Combining (CT ), { ^J% and ( ^TFD , we find that 

Zq = Z'^ J duj Z±(uj, uj') , Tq = T'^ + J dcj T±{uj, uj') . 
With our explicit expressions in ( ^.13| ), we obtain 



(3.17) 



(3.18) 
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To see in detail how this works, we calculate the ultraviolet divergences of the matrix 
elements {O\0±^'^'^'^'^\uj') to one-loop order, leaving the integration over until the end. 
Since now we are dealing with local operators, the only one-particle irreducible diagram is 
a vertex diagram analogous to that in Fig. Ha. Its contribution is: 



A;, 1 



(27r)2-2^ (fc2 ^ ^2 _ ,o)[A:2 + _ ^u') - zO] 
+ -Uu.^ IP^-^ 1, (3.20) 

where in the numerator of the second integral we have to substitute k- = k'j_/uj. This 
expression shows the origin of the two types of ultraviolet divergences. The k± integral in 
the first term is convergent, but there is a logarithmic divergence in the remaining integral 
over uj. In the second term, the u integral is cut off in the ultraviolet; however, the k± 
integral diverges for the minus projection. Keeping only the poles in 1/e, we obtain 

C^^7±«|±^"^/dc^^-'-'^ + (lTl)^ + ...| =C7^^7±«- (3.21) 

Adding to this the contributions from the wave-function renormalization of the external 
lines, we recover the expressions for Zq and Fq given in ( |3.14| ) and ( p.l5| ). The contribution 



in (|3.21|) arising from the logarithmic divergence of the u integral is removed by the factor 



Z'^ defined in ( |3.17| ), whereas all other contributions are removed by the renormalization of 
the operators 0±{u). 



IV. QCD SUM RULES 

The first moments of the wave functions ip±{uj) are known exactly from the equations 
of motion, as shown in ( p. 191) . However, they only set the scale of u. In order to obtain 
information about the shape of the wave functions, we need to consider some of the higher 
moments. In this section, we use QCD sum rules to investigate the parameters X% and A^, 
which according to ( |2.24D determine the second moments of ip±{u!). 



Let us consider the operators 

Oe = Q*id ■Eq = v.v'' Q^G^^'a^^ q , 

Oh = -Q*a ■Hq= - v,v") Q^G'^-'a^^ q , (4.1) 

whose matrix elements give /A|; and /A^, respectively. To obtain the sum rule for (the 
sum rule for Af^ is obtained in a similar way), we investigate the correlator 

{0\T{OE{x),qUl + ij)Q{0)}\0) = eiv x) 6{x±)Ue{v ■ x) , (4.2) 
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(c) (d) 

FIG. 4. Non- vanishing diagrams for the correlators II^; and Hh- The higher-dimensional cur- 
rent operators Oe and Oh are represented by a gray circle; the interpolating current is represented 
by a white circle. 

where x± contains the components of x orthogonal to v. The operator g |(1 -|- has the 
quantum numbers of the ground-state meson. Next, we analytically continue IIe{v ■ x) to 
V ■ X = —iT.f\ The correlator can be expressed via its spectral density Pe{i^)'- 

oo 

ns(r) = jdepE{e)e-'\ (4.3) 



which contains the contribution p\\5{e — A) from the ground-state meson, as well as 
contributions from excited states, which we will refer to as the continuum. Hence, the 
phenomenological expression for the correlator is 

^e{t) = fXle-^^ + nr (r) . (4.4) 

For sufficiently small values of r, the correlator can be calculated in QCD using the operator 
product expansion The theoretical spectral density, p%{e), contains perturbative as well 
as non-perturbative contributions, where the latter are proportional to vacuum condensates 
of local, gauge-invariant operators. The continuum contributions in ( |4.4| ) are usually mod- 
eled by the theoretical spectral density above a threshold Ec, which is called the continuum 
threshold. Equating the two expressions for IIe{t) obtained in this way, we derive the sum 
rule 

f\le-~^- = ldep'^{e)e-^\ (4.5) 



A similar sum rule holds for the product /^A|^. 



The procedure of analytic continuation of the coordinate-space correlator to imaginary time is 
equivalent to the Borel transformation of the corresponding momentum-space correlator. Thus, 
our sum rules coincide with the usual Laplace sum rules. 
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The leading contributions in the operator product expansion of the correlators are shown 
in Fig. ^. We include the leading perturbative contribution, as well as the contributions of the 
quark condensate {qq) (dimension (i = 3), the gluon condensate (G^) = {Ga^uGa^) (c^ = 4), 
and the mixed quark-gluon condensate {qa^^G^^q) = ml{qq) {d = 5). In the calculation of 
the non-perturbative contributions, it is convenient to use the fixed-point gauge x^A^{x) = 0; 



then the heavy quark does not interact with gluons. 
obtain the sum rules: 



where 



a. 



After a straightforward calculation, we 



16 ' 



ml{qq) 
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(4.6) 



5n{x) = 6{x) 



n ^in 



m=0 



ml 



(4.7) 



The fact that the sum rule for A|; does not contain contributions from the quark and gluon 
condensates is a consequence of the fact that, in the fixed-point gauge, the light quark 
interacts only with the magnetic components of the gluon field |53|. In Appendix B, we also 



list the contributions of higher-dimensional quark-gluon condensates {d = 7) to the sum 
rules. For the range of r values considered below, we expect that the contributions of these 
condensates are very small. Since their values are moreover unknown, they will be neglected 
in the numerical analysis. 

It is convenient to divide the sum rules in (|4.6|) by the sum rule for the meson decay 



constant derived in 42,18,20 



-At 



27r2 
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(4.8) 



This procedure leads to expressions for and Xjj as a function of r, the continuum threshold 



and the vacuum condensates. For our analysis we use the standard values ||5^ 



{qq) = -(0.23 GeV)^ , aAG^) = 0.04 GeV^ 



rrin 



0.8 GeV^ 



(4.9) 



as well as as = 0.4. The value of the continuum threshold extracted from the analysis of the 
sum rule ( [4.8| ) is Ec = 1.00 ± 0.15 GeV p!8|j20[| . Moreover, stability of this sum rule requires 



that the parameter r be in the range 0.3 GeV < 1/r < 0.5 GeV, which is called the stability 
window. Our numerical results for and as a function of 1/r are shown in Fig. ^ We 
observe a sizeable dependence of the results on the parameter r, as it is not unexpected for 
the matrix elements of higher- dimensional local operators such as and Oh- Nevertheless, 
taking an average over the stability window determined from the sum rule ( [4.8|) , we obtain 
as a rough estimate: 



X2 



(0.11 ±0.06) GeV^ 



X2 



(0.18 ± 0.07) GeV^ 



(4.10) 
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FIG. 5. Sum-rule results for (upper plot) and (lower plot) as a function of 1/r, for three 
values of the continuum threshold £c- The stability window lies in between the dashed lines. 



In units of A ^ 0.55 GeV this implies Al/A^ = 0.36 ±0.20 and A^/A^ = 0.60 ±0.23. 

According to ( |2.24] ), the parameters X% and X% determined, together with A^, the second 
moments of the meson wave functions. In order to learn more about the shape of the wave 
functions, we shall now go a step further and construct sum rules for the functions (p±{uj) 
themselves. To this end, we start from the correlator^ 



( |T{0±(t), g i(l ± ^)Q(-x)}| ) = 7± e{v ■ x) S{x^) U±{v ■ x, t) 



(4.11) 



where 0±(t) has been defined in ( |2.7D , and analytically continue Ii±{v ■ x,t) to v ■ x = —ir. 
Using the light-quark propagator in coordinate space , it is straightforward to obtain the 
contribution of the diagrams in Fig. ^,b. The result is: 



n«(r,t) 



2+2' 



2'kH{t + 2ity 



1 - 



^ Expanding this correlator in powers of t, one recovers, according to (2.13) and ( 2. 151 ), the sum 
rules for the moments f'^{uj'^)±- To lowest order, this gives the sum rule for the decay constant 

The sum rule for /"^{lo) 



f considered in [42,15 2C|. The sum rule for f'^{uj)^ has been considered in ||5j]; however, the 
mixed-condensate contribution had the wrong sign. Note that, because of (pl9| ), the sum rules 
for f'^{u>)± can be obtained by taking derivatives with respect to r in the sum rule for Taking 
linear combinations of the sum rules for f'^{u>'^)±, we recover the sum rules for A|; and X^. 
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FIG. 6. Leading diagrams for the correlators of the bifocal operators 0±{t) (gray circles) with 
the local current qQ (white circle). 



n«(r,t) 



1 - 



(4.12) 



The contribution with the cut light-quark line shown in Fig. ^ involves the trilocal, non- 



collinear quark condensate This object is discussed in detail in Appendix B. Keeping 
for simplicity only the leading term in the operator product expansion, we obtain 



1 



{qq) / dz//s(z/)e 



-l/T{T+2it) 



(4.13) 



where /s'(z^) describes the distribution of quarks with virtuality u in the vacuum. Unfor- 
tunately, little is known about the shape of this function. A simple ansatz is discussed in 
Appendix B. 

Next we perform the Fourier transform of the correlator in ( [4.11 ) with respect to t, which 
leads to 



|T{0±(cu), g i(l + ^)Q(-x)}| ) = 7± 



1 + ^ 



{v ■ x) S{x±) Il±{v ■ X,U!) . 



The leading perturbative terms in ( [1.12| ) give 



n«(r,o;) 



u e 



-u)t/2 



-ujt/2 



(4.14) 



(4.15) 



The contributions proportional to the gluon condensate lead to singular behaviour of the 
form S{uj) and 6'{uj). These terms would acquire a finite width if the non-locality of the gluon 



condensate would be taken into account [0]. Below, we shall neglect the gluon condensate. 
The Fourier transform of the quark-condensate contribution in ( [4. 131) is given by 



n(2) 



{qq) 

8r 



fs 



2tJ 



r/2 



(4.16) 



i.e. it is directly determined by the virtuality distribution of quarks in the vacuum. 

In order to perform the continuum subtraction for the perturbative contributions, we 
calculate the inverse Laplace transforms of the expressions in (|4.15|) , which give the corre- 
sponding spectral densities. We finally arrive at the sum rules 



/V-(^)e-^^ 



SttV 

4^2^2 



6l [(^c-f)T 



(qq) 7 , 



ujt/2 
-ujt/2 



(4.17) 
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FIG. 7. Model wave functions ip±{uj) defined in ( 4.1 



The contributions from the non-local quark condensate fall off quickly for both u; ^ 
and u oo. The first statement follows from a general property of the function fs (see 
Appendix B). However, the precise functional form of these contributions is unknown. The 
leading perturbative contributions vanish for uj > 2ec- They suggest the model shapes 



UJ 
^0 



-uj/u)o 



(4.18) 



which do indeed exhibit the correct behaviour for c<j — > [cf. (|2.33|) ]. These functions are 



shown in Fig. |^. The parameter uq = |A is fixed by (|2.19|) . In this simple model, we obtain 



for the second moments {uj'^)+ = 3{uj'^)- = |A^. This corresponds to 
does not contradict our sum-rule estimates obtained earlier in this section. 



X2 



|A^, which 



V. ASYMPTOTICS OF FORM FACTORS 

We shall now use the results obtained in the previous sections to analyse, in a model- 
independent way, the asymptotic behaviour at large recoil of the form factors describing 
the current matrix elements between two heavy mesons. The contribution of the quark- 
antiquark wave functions to the Isgur-Wise form factor is depicted by the diagrams in 
Fig. |^a,b. To deal with two heavy mesons moving at different velocities v and v', it is 
convenient to choose the Breit frame, in which the two mesons move in opposite directions 
with rapidities ±-(9/2, so that 

i;'' = (cosh ^,0,0, sinh^), t;'^ = (cosh ^, 0, 0, - sinh ^) , (5.1) 

and V ■ v' = coshi?. In terms of the light-cone vectors n±, we have 

w/^ + t;'/^ = cosh-« + n^), -v'^ = -smh-{n'^-n1) . (5.2) 

It follows that f+ = e'^/^ and V- = e~^^'^, but v'_ = e^l"^ and v'j^ = e~^^'^. Similarly, in 
the large-recoil limit {'& ^ 0) the light quark in the initial meson has large p+ = ue^^^ 
and small p_, whereas the light quark in the final meson has large p'_ = Lo'e^l"^ and small 
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FIG. 8. Hard-gluon exchange contributions to heavy-meson form factors. The external current 
is represented by the wave hne attached to the gray circle. The white square in (d) represents an 
insertion of l/mg-suppressed operators from the effective Lagrangian of the HQET. 



Thus, the roles of the plus/minus directions for the final meson are opposite to those 
for the initial meson. The virtuality of the gluon and of the heavy quark are both large: 
k"^ = {p' — pY — —ujuj'e^, and v ■ kq c:^ v ■ p' (Fig- IK) or t>' • A;q ~ t>' • p ~ 

(Fig. 1^), respectively. 

For large the contribution of the diagram in Fig. §a to the matrix element in ( p.. 41) , 
which defines the Isgur-Wise function, is 



2™. — / e 



LULU ■ 



-ui 



X 



V9+(cj)7_e''/^ + v2-(u;)7+e~''''^l u{v) 



(5.3) 



where f = ^ {-f.e^/^ + -f+e-'^/^) . The factor arises from the normalization of the 



colour wave functions of the meson states, which is such that the matrix element in (|LT 
normalized to the decay constant /. Using the relations in 



IS 



together with ( p.2| ), we find 
that the leading contribution contains the product ip+i^u') (p-^uj), i.e. the subleading-twist 
wave function is taken on the side where the gluon exchange occurs. Adding the contribution 
of the diagram in Fig. Hb, we arrive at 



^(cosh??) = 16na. 



(5.4) 



Based on our assumptions about the behaviour of the wave functions for — > 
[cf. ( |2.33|) ], we expect that both (0;"^)+ and (a;~^)_ diverge logarithmically at low u. This 
divergence is cut off by the transverse momenta and virtualities of the light quarks in the 
mesons, similar to the case of the vr-p form factor in QCD ||29[. This infrared sensitivity 
results in an additional enhancement of the form factor, as can be seen by replacing by 
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g2 _ ^YiQ gluon propagator (where A is of the order of a typical hadronic scale) , which 
leads to the replacement: 

Unfortunately, however, the subleading terms of order d cannot be calculated without know- 
ing details of the infrared cutoff. 

Consider now the form factors E,i and Xi defined in ( |1.5| )-( p^ ), which appear at order 
1 / niQ in the heavy-quark expansion. Since their contributions to the meson form factors are 
suppressed by a power of A/mg with respect to the contribution of the Isgur-Wise function, 
they can only become important if they have a slower fall-off at large d. Therefore, it is 
sufficient to retain only those contributions where the leading-twist wave function appears on 
both sides of the diagrams. Then the meson helicity is conserved, and the gluon polarization 
is orthogonal to the v-v' plane. Let us first focus on the functions ^i. The only way to get the 
leading-twist wave function on both sides of the diagram is to attach the gluon to the current 
operator, as shown by the diagram in Fig. |^c. Simplifying the resulting spinor product using 



]8|) and ( |5.2| ), we obtain 



4™, ^ p u{v')Yu{v) f ^4^) Mu^') , (5.6) 



Nc J UJUJ' 

which means that only the function ,^3 receives a leading-twist contribution. It is given by 

Ucosh{}) = 4nas^ f {uj-%e~^ . (5.7) 

Here the integral is infrared convergent, as in the case of the pion from factor in QCD 
25|^T1. Next consider the functions Xi- The conservation of the meson helicity implies that 



there must be an odd number of 7 matrices in the matrix element. Indeed, only the function 
X2 in ( |1.9| ) receives a leading contribution. Calculating the diagram in Fig. |^d, we find 

- 8na, ^ f e-'' f ^^{uj) ^^{uj') u{v'){rkl - i^k';)u{v) , (5.8) 

1\ Q J UJUJ 

where kg = uj'v' — ujv is the gluon momentum. Since, by definition, the indices /i and v are 
restricted to the subspace orthogonal to f , only the first term in kg has to be kept. Taking 
into account the definition of X2 in ( |1-8|) and (|1.9| ), we obtain 



X2(cosh^9) = -^3(cosh79) e"'^ . (5.9) 

Note that in the expressions for the meson form factors X2 is multiplied by cosh d > so 
that its contributions are of the same order as the contributions of ^3. 

We can compare our asymptotic results for the leading and subleading Isgur-Wise func- 
tions in (p.4|), ( |5.7| ) and ( |5.9| ) with the large- limit of the two-loop QCD sum-rule expressions 



for these functions, which have been obtained in pl] , p4| . We find that the results of the sum- 
rule calculations do indeed reproduce the correct asymptotic behaviour; in particular, the 
relation (|5.9|) between X2 and ^3 is satisfied. Moreover, the sum rules allow us to determine 



the normalization factors appearing in the expressions for ^ in (|5.4|) and for ^3 in (|5.7|). This 
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is explained in detail in Appendix C. For later convenience, we also present the expressions 
obtained using the model wave functions (|4.18|) . They are: 

P ln^(v ■ v') P 1 
i{vv')^^7,a,L-^ i,{vv')^27:a,L^ -. (5.10) 

Let us discuss the applicability regions for these asymptotic results. QCD sum rules suggest 
that there are "soft" contributions to the Isgur-Wise function which fall of like l/(f ■ v'Y 
p^7|^0[. If this is correct, the asymptotic behaviour given by ( |5.4| ) and ( p. 51 ) would dominate 
only if ln(f ■ v') S> 1 and as\v?{v ■ v') S> 1. If the "soft" contributions vanish faster than 
l/{v ■ v'Y, the second requirement is removed. The fact that ln(f ■ v') is, in most practical 
applications, not a large parameter implies that the asymptotic result for the Isgur-Wise 
function may be considered as a rough estimate only. To reach the asymptotic regime would 
require v ■ v' = O(IOO). QCD sum rules also suggest that the "soft" contributions to ^3 
fall off like l/{v ■ v'Y p2|- p^] , meaning that the leading hard contribution given in ( p.7[ ) is 
enhanced by a power of v ■ v' . As a consequence, our predictions for the functions ^3 and 
X2 are much more accurate than for the Isgur-Wise function. The asymptotic behaviour 
should set in when f ■ f ' ^ 1, which requires v ■ v' = 0(10). 

An important aspect of physics is still missing from our discussion of the asymptotic 
behaviour of meson form factors at large recoil. Since the quarks receive a large accel- 
eration during the transition process, they emit gluon bremsstrahlung, which leads to an 
additional damping of the transition amplitudes (Sudakov form factor). Because the mesons 
are colourless, the double logarithms of the type [ug \v?{v ■ f ')]" cancel in the expressions for 
the meson form factors; however, single logarithms in v ■ v' remain, which are enhanced by 
logarithms of the heavy-quark mass. They arise from the emission of gluon bremsstrahlung 
with energies in the range n < Eg < mq (/z -C mg, see below). Thus, in perturbation theory 
there are large double-logarithmic contributions of the type [a^ ln(mQ//i) ln(u • v')Y to the 
form factors. The situation is similar to the case of the contributions to the pion form factor 
coming from the region a; — > 0, where almost all of the pion momentum is carried by one 
quark |5^ . 



Because of the explicit dependence on the heavy-quark mass, these large logarithms are 
not contained in the form factors of the HQET, which are renormalized at a scale /i -C 
However, they appear when we relate the form factors of the HQET to physical meson 
form factors using a perturbative matching procedure. In this relation, there appear short- 
distance coefficient functions Cn{fnq/ fi,v ■ v'), which can be calculated in renormalizat ion- 
group improved perturbation theory (see 0] for a review). In the case of transitions between 
two heavy hadrons, the reparametrization invariance of the HQET ensures that the 



coefficients multiplying the subleading functions X2 and ^3 are the same as the coefficients 
multiplying the leading-order Isgur-Wise function C, Indeed, in leading logarithmic 

approximation (which is sufficient to control the large logarithms mentioned above), all 
HQET form factors are multiplied by a universal coefficient 



^In practice, the scale fj, should be chosen such that there are no large logarithms contained in 
the form factors of the HQET, but yet large enough for perturbation theory to be valid. A typical 
choice is /i ~ 1 GeV. 
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Cimg/n^v v') 



a{v-v') 



(5.11) 



where {v ■ v' = cosh -d) p5|j60|j6T[] 



a{vv') = =^{^coth'd-l) ^ 



For large recoil, we find that 



^-l + 0(e- 



(5.12) 



C 



2v ■ v' 



where 



V 



20 F asifi) 



(3o asimq) 



(5.13) 



(5.14) 



This expression sums the large Sudakov logarithms correctly to all orders in perturbation 
theory. The effect of this bremsstrahlung correction is an additional power-like suppression 
of the physical meson form factors for large values ofv-v'. Using fi ^ 1 GeV, we find that for 
S-meson decays the power rj is given by 77 ~ 0.2-0.3, i.e. the overall effect of bremsstrahlung 
emission is rather small. 

Using the results obtained in this section, it is straightforward to derive the asymptotic 
behaviour of all form factors describing current-induced transitions between any two pseu- 
doscalar or vector mesons containing a heavy quark. The relevant formulae, which relate 
the meson form factors to the Isgur-Wise functions, can be found, e.g., in P JT6|J51|] . Here 
we restrict ourselves to the results obtained for the matrix elements of the vector current 

= b-f^b between S-meson states. We find^ 



{B{v')\V^'\B{v) 
{B*{e',v')\V^'\B{v) 
{Bl{v')\V^\Bl{v) 
{BUe',v')\V'^\BUe,v) 



(B* 



[e',v')\V'\Bl{v: 



h+ {v + vY , 

hL {v + v'r, 

-/iT e ■ e* {v + vY , 

sinh^^/iTL e'*^', 



(5.15) 



where L and T refer to longitudinal (i.e. in the v-v' plane) and transverse (i.e. orthogonal 
to that plane) polarization states. We find that, asymptotically. 



C 



nib 
1 



cosh ^ X2 



rub 



e- — (e3 + 2cosh^X2) 
nib 



To obtain the conventional relativistic normalization of meson states, the right-hand sides in 
these equations have to be multiplied by 
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hL = cU + — cosh^^xa) =CU es) 



rrib 



C^. (5.16) 



The most striking feature of these resuhs is the fact that the form factor for two longitu- 
dinally polarized B* mesons, which is positive for cosh?? <C rrib/A, becomes negative for 
cosh-t? ^ rrih/A, and hence has a zero at some intermediate values cosh-;? ~ mf,/A. Since 
for spacelike (negative) values of g^, corresponding to f ■ f' > 1, all form factors are real, 
the existence of this zero is an exact statement not affected by subleading corrections. We 
should stress that this observation is not specific to heavy-light mesons. The form factor 
of, say, longitudinally polarized p mesons also has a zero at some negative value of p9| . 



For timelike (positive) values of g^, on the other hand, it is the form factor of pseudoscalar 
B mesons which has a zero, and this zero is situated inside the physical region of the pro- 
duction of BB pairs in e~^e~ collisions. Strictly speaking, because form factors at timelike 
values of g^, corresponding to f ■ f' < — 1, are complex, this zero is not absolutely exact. 
However, in our approximation the imaginary part is negligible.P^ 

The model-independent results obtained in this section can be checked in the simple 
model where a meson is built out of two heavy quarks with masses m and /i such that 
m ^ /i ^ A. This is discussed in detail in Appendix D. The same model has been considered 
by Brodsky and Ji who observed for the first time the zero of the pseudoscalar form 



factor in the physical region of large positive g^. However, their claim that the form factor 
of longitudinally polarized vector mesons would have the same behaviour is incorrect. Our 
analysis shows that this form factor has a zero in the region of spacelike momentum transfer, 
i.e. for large negative g^. 

VI. APPLICATIONS 

We finally apply our results to calculate the cross section for the reaction e"'"e~ ij(*)^(*) 
in the region s ^ 4m^. Recall from the introduction that the part of the electromagnetic 
current that couples to light quarks does not give a leading contribution in the asymptotic 
regime. Hence, it is justified to use the relations in ( |5.16| ) for the relevant form factors. As 
usual, we define 

n^.^if^^ (6.1) 



Using crossing symmetry and the matrix elements given in ( |5.15| ), we obtain for x 
^IniB > 1: 



^'^Note also that any mechanism that leads to a common phase factor of the form factors, such as 
final-state interactions, does not spoil the existence of the zero. 
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Ri 


= R 


R2 


= R 


R3 


= R 


i?4 


= R 



BB = \\C\' 



2 2 ^ 

z2 2 2 



5 



e 6 

rub 



j\Cm', (6.2) 



where the form factors are functions of f ■ f ' ~ — l^;^, and \C\ = (e/x^)^'' with r] given in 
( ^■14| ). Here Zf, = — | is the electric charge of the b quark. As an illustration, we show in 
Fig. 1^ our predictions for the various B^*^B^*^ production cross-sections as a function of x 
obtained using the model wave functions in (|4.18| ), for which [cf. ( p.lO|) ] 



e(a:)~ 487r — ^3(3;) ^ -47r — - . (6.3) 



For simplicity, we have neglected bremsstrahlung effects setting C = 1 in (|6.2|). Because 
of the kinematic enhancement factor x^, the ratio -Ri generally dominates at large x. This 
means that mostly BB* and B^B^ pairs are produced, with i? ~ and angular distribu- 
tion 1 -|- cos^ 6 ||62|,^ . More interesting from the point of view of the present work are the 



three other ratios. The BB production cross section vanishes at some value Xq ~ ymb/A, 
since ^ ~ In^x/x^ and .^3 ~ — A/x^. In our simple model, we find that: 



xo ~ y ^'^^o ~ 25. (6.4) 

We stress again that the accuracy of this prediction is not high, because of our poor knowl- 
edge of the asymptotic behaviour of the Isgur-Wise function. For 1 ^ x <^ Xq, the ratios 
R2, R3 and R4 are all of the same magnitude and scale like 1/x^. This situation had been 
studied previously in the context of the HQET |]62| , |63| . For x ^ xq, however, another 



pattern sets in. Then the contribution of ^3 to the ratios R2 and Rs dominates over the 
contribution from the Isgur-Wise function ^, so that R2 and R3 scale like 1/x^ and domi- 
nate over R4. In principle, at very large x the ratios R2 and R^ should even dominate over 
Ri, which scales like 1/x^. However, because of the double-logarithmic enhancement of the 
Isgur-Wise function this would require enormous values x > xi, where in our model xi is 
given by 

x-^ ^ ^ In^xi ^ 3500 . (6.5) 

In the ultra-asymptotic region x 3> xi, mostly BB and B^B^ pairs would be produced, with 
i? ~ and angular distribution sin^ 

These qualitative features, which are independent of the particular choice adopted for 
the meson wave functions, are clearly exhibited in Fig. ^. Unfortunately, however, the 
cross sections for B^*^B^*^ production at large x are so small that they will most likely be 
irrelevant to experiments. The situation is somewhat more favourable in the case of the pair 
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FIG. 9. Cross-section ratios Ri for BB pair production in e~^e collisions at large x. We use 
the form factors in 
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FIG. 10. Cross-section ratios i?j for DD pair production in e^e collisions at large x. 



production of charm mesons. We can apply our results to this case by performing obvious 
substitutions (m?, rric, ^ Zc = |, etc.) in the above formulae. We then find that xq ~ 8 
in the case of charm pair production, corresponding to moderate energies of order 15 GeV. 
The second turnover point is, however, still too high (xi ~ 650) to be of any interest. The 
resulting cross sections are shown in Fig. 

In summary, we have applied methods developed for hard exclusive QCD processes to 
calculate the asymptotic behaviour of heavy-meson form factors at large recoil. We find that 
this behaviour is determined by the leading- and subleading-twist meson wave functions. For 
1 \v ■ v'\ <^ ttiq/A, the form factors are dominated by the Isgur-Wise function ^, which 
is determined by the interference between the wave functions of leading and subleading 
twist. At It" ■ f'l ^ mg/A, they are dominated by the two functions ^3 and X2 arising at 
order l/mq in the heavy-quark expansion, which are determined by the leading- twist wave 
function alone. The sum of these contributions describes the form factors in the whole 
region \v ■ v'\ ^ 1. Central objects of our study are the meson wave functions (f±{uj), which 
are defined in terms of the Fourier transforms of the matrix element of bilocal operators 
on the light-cone. We have derived the (Brodsky-Lepage) evolution equations obeyed by 
these wave functions, and we have investigated the properties of the wave functions (such 
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as their moments) using QCD sum rules. Finally, we have discussed as an application the 
implications of our results for the production of heavy-meson pairs in e^e~ collisions. 
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APPENDIX A: COVARIANT TRACE FORMALISM 

The most convenient way to calculate the matrix elements of operators between the 
physical pseudoscalar and vector (qQ) meson states (rather than the spin-^ mesons obtained 
when the heavy-quark spin is switched off) is provided by the covariant tensor formalism 
introduced in [|l5| . In the HQET, the spin wave function of the ground-state meson doublet 



is described by the 4x4 Dirac matrix 

_ ^ ~^ i ""^75 ) pseudoscalar meson M(t>), (Al) 
2 1 ^ ; vector meson M*(e, f), 

where v is the meson velocity, and e is the polarization vector of the vector meson {e-v = 0). 
The matrix /iA{v) simply contains the appropriate spin-parity projections of the spinor 
product Ug{v) vq{v) [Q. It satisfies 

^^M{v) = M{v)M{v) = -Miv)^. (A2) 

Operator matrix elements between meson states can be represented by traces over these 
wave functions. Consider first the matrix elements of heavy-light operators of the type 
QyT 0{iD) q, where F is an arbitrary Dirac matrix, and 0{iD) is a differential operator 
acting on the light-quark field, between a meson state and the vacuum. Their representation 
is 

( \Q,T OiiD) q\M{v)) = Ti{Oiv)Miv) F} , (A3) 

where 0{v) is the most general matrix with the same transformation properties (under the 
Lorentz group and heavy-quark symmetry) as the operator O. In the spinor formalism 
adopted in our paper, the same matrix element would read 

{Q\QlO{iD)q\M{v)) = Oiv)uiv) (A4) 

with the same matrix 0{v). The covariant decomposition of this matrix determines the 
number of reduced matrix elements (generalized Isgur-Wise form factors) that appear in 
the heavy-quark expansion. As an example, we give the expressions in the trace formalism 
which correspond to the definitions in (|1 . 1| ) , ( p. 181) , and (|2.20| ): 
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{0\Q,T q\M{v)) = fTi{M{v)T} , 
{0\Q^riD''q\M{v)) = |/ATr{(4f'' -7^)A^(t;)r} 
'0\Q^riD^'iD''q\M{v)) = fTT{Q^"'{v)M{v)r} , 



(A5) 



where B'^'^(f) is given in ( p.21| ). 

In the end of Sec. 0, we need the generahzation of ( |2.3| ) in the trace formahsm. It reads: 



|Q,(0)r EiO, z) q{z)\M{v)) = f Tr{ + i - Miv) v] . (A6) 



Evaluating the trace for various choices of F, we recover the results given in ( |2.26D -( p^.31D . 

The trace formalism is readily extended to more complicated cases, such as transition 
matrix elements between two meson states. For instance, the expressions corresponding to 
the definitions in (|1.4|) and ( |1.5|) read: 



{M{v')\Q,TQ,,\M{v)) = ■ v') TT{Miv)T Miv')} , 
{M{v') I {iD''^Q,)TQ,,\M{v)) = Trieiv, v') M (v) F M{v')} 



(A7) 



where the covariant decomposition of ^'^{v, v') is given in ( [I.6| ). Note that as a consequence of 
the fact that in this paper we work with (qQ) rather than (Qq) mesons, the trace formalism is 
slightly different from the one usual employed in the literature P,p!6|j51| . Crossing symmetry 



implies that the form factors for (qQ) mesons are related to those for (Qq) mesons by 
Hermitean conjugation followed by the substitutions v —v and v' —>■ —v'. We have 
defined the invariant functions in the HQET (^, ^3, Xi etc.) in such a way that the resulting 
expressions for the physical matrix elements look the same as in the conventional formalism. 



APPENDIX B: NON-LOCAL CONDENSATES 

The contributions of higher-order non-perturbative corrections to the sum rules for A|; 
and X'jj in (|4.6| ) can be included by introducing two functions, f^^\x'^) and f^'^\x'^), which 
parametrize the following non-local condensates |]55| : 



2 ^ (g(0)E(0,x)a^.G^'-(x)g(x)) ^ Qi - Q2 - 2Qs ^ 
^ ^ {qa^^G^^-^q) ^ ml{qq) 9> ^ ' ' ' ' 

(2). 2^ _ 4x"a;^ (g(0)E(0,x)a^^G^^(x)g(x)) 2Qi - - SQa , . 

^ ' x2 (ga^.GA'-g) + mliqq) 12 + ' ' ' ' ^ ^ 

The quantities Qi form a basis of dimension-7 quark-gluon condensates and are defined as 

Qi = (qG^.G^'q) , 
Q2 = i{qG^,G^''l^q) , 
Qz = i{qo,uG^'^G\q), 

Q,= {qa^,iD'^D^G'''')q). (B2) 
If these corrections are included, the sum rules ([4.6|) are modified in the following way: 
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ml{qq) 



16 



ml{qq) 



IGttt 



16 



2;(l)(_^2)_^(2), 



(B3) 



Next we give some details of the calculation of the quark-condensate contributions to the 
QCD sum rules for the meson wave functions. The contribution with the cut light-quark 
line shown in Fig. ^ involves a trilocal object, the non-coUinear quark condensate [^] : 



{qp{y)E{y, 0)E(0,x)g,(x)) 



fsix, y) + ^ ^] fAx, y)- - [it fvix, y)-ij fviy, x)] 



(B4) 



J a/3 



Neglecting the function fv{x,y), whose operator product expansion contains even-dimen- 
sional quark condensates with d > Q and whose contribution to heavy-meson sum rules is 
neghgible |T8|,|40|-p| , we obtain 



n±V -x^t) = -i {qq) fs{z, -x) ± ^mlv ■ xt friz, -x) 



(B5) 



where z"^ = 0, x^ = {v ■ xY, and z ■ x = v ■ xt. 

The non-coUinear condensate in (|B4|) can be expanded in x at fixed y. One finds 



fs{x,y) = fs[{y-xf] + 
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4 (xy)^ — x^y"^ 



r 



rrin 



+ 0{x^ 



xy 



fTi^^y) = —2\fsiy')--f 



fW) + ^y'fsiy')-^f^'\y') 



+ 0(x' 



(B6) 



where the function fsix"^) parametrizes the bilocal quark condensate and is given by |]43|j55 



2^ {mEiO,x)q{x)) ^ , mix' , 6Q^ - 3Q2 - 6Q3 + 2Q, x^ ^ ^, 

JS[X ) = 7— T = i H — \ 7— T 7-377 + U[X ) . [ti/) 



{qq) 



16 



{qq) 



1152 



A convenient representation of the bilocal quark condensate is 

fs(x')= Jdufs{u)e^^\ 

where 

'diyfsiiy) = l 



rrin 



(B8) 



(B9) 



and so on. The function fs{i^) can be interpreted as the distribution quarks with virtuality 
u in the QCD vacuum. The local expansion in (P7|) corresponds to the expansion 



rrin 



(BIO) 
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Because the factor E{0, x) in (p7| ) can be interpreted as the heavy-quark propagator, the 
asymptotic behaviour of fsix"^) at large — is 



fs{x') e-^^-^\ (Bll) 

This fixes the behaviour of fs{i^) for i> ^ 0. A simple ansatz for the distribution function, 
which satisfies this constraint, was proposed in ||64|| : 



4z/ 

where N and a are fixed by the conditions ( p9|) 



fs{u) = Nexp(-^-au] , (B12) 



APPENDIX C: FORM-FACTOR ASYMPTOTICS FROM QCD SUM RULES 

It is instructive to compare our asymptotic results for the leading and subleading Isgur- 
Wise functions with the large-t? limit of the two-loop (order-a^) QCD sum-rule expressions 
for these functions, which have been obtained in pT|J2^ . For very large recoil, the three- 
point correlators considered in these sum rules factorize into the convolution of the two-point 
correlators ([4.14|) with hard-scattering amplitudes. In this limit, only diagrams with a gluon 



exchange between a heavy quark and the light quark remain. We find that the results of the 
sum-rule calculations do indeed reproduce the correct asymptotic behaviour; in particular, 
the relation ( |5.9| ) between X2 and ^3 is satisfied. For the normalization factors appearing in 
the expressions for ^ in ( |5.4|) and for ^3 in ( ^.7|) , we obtain from 



Ar2 

/'(--^)l = ^^'^^3(.cr). (CI) 

We have retained the leading perturbative contributions only, since the relevant non-local 
condensates have not yet been calculated to order ag- Note, in particular, that the leading 
quark-condensate contribution to the sum rule for the Isgur-Wise function is constant and 
seems to dominate for large recoil. However, once the non-locality of the quark condensate is 
taken into account, one finds that this contribution actually vanishes quickly at large recoil 
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It is straightforward reproduce the expressions in ( PT| ) starting from the sum-rule results 



for the wave functions ^±{00) obtained in Sec. |^ The sum rule for the product (p^{uj)(p±{oii') 
at equal Borel parameters has the form 



where p±{uj,e) are the spectral densities of the correlators ( 4.14 ), and the integral is taken 
over the complement of the continuum region. The precise form of the result will depend on 
the particular way in which the continuum subtraction is performed. The spectral densities 
for three-point correlators depend on two variables, e and e'; see Fig. |lT|. The "square model" 
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FIG. 11. Square model (solid line) and triangle model (dashed line) for the continuum subtrac- 
tion in QCD sum rules for three-point correlation functions. 

of the continuum subtraction amounts to cutting off the integrals over these variables at 
the threshold Sc- This leads to an exact factorization of the integrals in ( |C2|) . Using the 
appropriate products of the sum rules for the wave functions given in ( ^.17| ), and retaining 
the leading perturbative contributions only, we then obtain: 



/V+(o)^-(o) 



(C3) 



On the other hand, it is well known that for small recoil the square model of the continuum 
subtraction is inconsistent, as it leads to an unphysical infinite slope of the Isgur-Wise 
function at f ■ f ' = 1 |18|. This deficiency is removed by using the "triangle model", where 
< e + < 2eci while the difference e — e' is unconstrained ||T8| , p!9[| . The triangle model 
was adopted in the calculations in pl| , |2^ . If we use it to evaluate (|C^), we indeed recover 
(pll). This is a strong check of both, the present approach and the two-loop calculations 
performed in However, in this case the resulting sum rule for the product of the 

wave functions is no longer exactly factorizable, meaning that the triangle model is not fully 
consistent at large recoil, and so the square model is preferrable. Comparing ( PT| ) with (|C^), 
we observe that both results agree in the limit Sc ^ oo, when the choice of the continuum 
model becomes irrelevant (5„ — >■ 1). 



The sum rule for the Isgur-Wise function |T8| has a built-in model of the infrared cutoff, 
and therefore it allows us to estimate the subleading 0((d) term to the Isgur-Wise function 
in (15.51). The result is 



where 



e(cosh d) ~ e-^^ [d"^ + [4L(£er) - 5] ?9 + . . . } , 
/ dxx \n2x 



LiXr 



J dx x'^e~^ 





In the relevant region of values of EcT ~ 2.5, we find that (4L — 5) 
the subleading term is negative and has a coefficient of order unity. 



(C4) 
(C5) 

-0.7, meaning that 
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APPENDIX D: STATIC QUARK MODEL 



Our model- independent results in {f>A\), (|5.7D, (|5.9|), and ( |5.15| ) can be checked in the 
simple model where a meson is composed of two heavy quarks with masses m and fi, such 
that m ^ yU ^ A [|^. In this case, (p+{uj) = (p-{uj) = 6{lj — /i), and {uj"')± = /x". Note that 
the formulae ( p.l9|) and ( |2.24|) for the lowest moments are based on the equation of motion 
for a massless quark, and are thus no longer applicable. 

Let us consider the pseudoscalar form factor in ( |5.15| ) . It is convenient to calculate it 
from the relation (cosht?-!- 1) = {B{v')\v^V'^\B{v)) . A simple evaluation of the diagrams 
in Fig. Ra,b gives 



2nas ITT ^ e Tr 



(Dl) 



where $ = 75(1 — and $' = —(1 — i)')'^^ are the spin structures arising for pseudoscalar 
mesons. The heavy-quark propagators in Fig. |]a,b are given by 



S 



m{l — i^) + l^f 
— m/ie"^ 



(D2) 



and S' = S{v ^ v'). If we retain the leading term proportional to m in the numerator 
of (|D2| ) , then the gluon is longitudinally polarized; the two diagrams contribute equally, and 
we recover the contribution of the function C, in ( p.4| ) to the form factor in ( |5.16| ). If, on 
the other hand, we retain the term with we lose a factor /i/m but gain a factor e"^ from 
the trace. Then only the first diagram contributes; the gluon is transversely polarized, and 
we recover the contribution of the function ^3 in (|5.7| ) to the form factor In this model, 
the form factor of pseudoscalar mesons has a zero at = 2m^//i. The calculation can be 
repeated for vector mesons using the spin structures $ = ^(1 — and $' = (1 — In 
particular, we find that the form factor of longitudinally polarized vector mesons has a zero 
at = —2m^/fi. Other form factors can be calculated in a similar way; our results agree 

with Eel. 
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